INTEGRALES 
MULTIPLES 


por 
MIGUEL DE UNAMUNO ADARRAGA 


CUADERNOS 
DEL INSTITUTO 
JUAN DE HERRERA 
DE LA ESCUELA DE 
ARQUITECTURA 
DE MADRID 


3-07-02 


CUADERNOS 
DEL INSTITUTO 
JUAN DE HERRERA 


0 VARIOS 
Il. ESTRUCTURAS 
CONSTRUCCIÓN 

FÍSICA Y MATEMÁTICAS 
TEORÍA 

GEOMETRÍA Y DIBUJO 
PROYECTOS 
URBANISMO 
RESTAURACIÓN 


0 30 (Un hh u NN 


NUEVA NUMERACIÓN 


3 Área 
7 Autor 


02 Ordinal de cuaderno (del autor) 


Integrales múltiples 

(O 2002 Miguel de Unamuno Adarraga 

Instituto Juan de Herrera. 

Escuela Técnica Superior de Arquitectura de Madrid. 
Gestión y portada: Laura Bejerano Iglesias 
CUADERNO 15.02 / 3-07-02 

ISBN: 84-9728-083-0 

Depósito Legal: M-35902-2003 


I - Integral doble 


Integral doble de una función acotada en un rectángulo 
Consideremos una función acotada 
fiR>R, 
siendo R =|[a,b]x[c,d] un rectángulo del plano real, y sendas particiones 
G=ÍY4 = 6, senor Y =$ 
de los dos intervalos [a,b] y [c,d] 
de R. Quedarán definidos m-n 


rectángulos elementales (en la 
figura m=7,n=6), 


KR, =[x, ,, x,]xly,.,. »,), 


en cada uno de los cuales /, como 
5 función acotada, tendrá dos 

6 

extremos, 


M,,=sup f(x y), m,=inf f(x, y). 
i,) 1] 
Formemos las sumas 


s(p.q)= 2 M,, A(R,,). s(p,q)= 2, a(R,,), 


donde A(R, ,) es el área de R,,. 


Evidentemente, 


s(p,9)<S(p, q). 


Y lo mismo que en una variable, se demuestra que, para cualesquiera particiones p, p” 
de. [a,b] y q, q' de [c,d], es 


s(p.q)<S(p',q'). 
Con lo cual, 
DE Ís (p.q)| Pp, q particiones de [a, b), [c, d ] resp.) 
es un conjunto inferiormente acotado (por cualquier s(p,q)), y 


o=4s(p, a)! p, q particiones de |a,b), [c, d ] resp.) 


es superiormente acotado (por cualquier Sí p,q) ); luego existen su ínfimo y su supremo, 
respectivamente. Llamamos entonces: 


- Integral superior de f en R: S, f =inf Z,e integral inferior de fen R: Í. f=supc0, 
siendo evidentemente 


=> 


PL 


Definición. Se dice que fes” integrable en R si y sólo si 


a este número real se le llama entonces integral de fen R, con las notaciones 


1,£. O |, Sad, O JA y)ardy, O bien ff, vaca, 


pudiéndose escribir también da en lugar de dx dy. 
Ejemplo 1 


Sea f(x,y)=k e R una función constante en un rectángulo R la b]x[c,d], y p y q 


dos particiones cualesquiera de [a,b] y [c,d]. Será M, ¡=m,,=k para todos los 


valores de ij, con lo cual S(p,q)=s(p,q)= eee fs = ¡E f es 
integrable y 


|, 7 =k(b=a)(d-c) 


(si k > 0, es el volumen del prisma recto de base R y altura k). 


Ejemplo 2 
Sea ahora la función f: R => R definida así: 


sl si (x,y)e 0”, 


O en otro caso; 


con toda facilidad comprobaríamos que, cualesquiera que sean p y q, serán S (p, q) = 


= (b-a)(d-c)= J.7, y s(p,q)=0= f,, : luego fno es integrable en R. 


La integral doble como límite 


Análogamente al caso de la integral simple, podemos definir la integral doble como 
límite, mediante el concepto de suma de Riemann. 


Dadas como antes dos particiones p y q de los intervalos [a,b] y [c,d], llamamos suma 


de Riemann correspondiente a las mismas a plp,q )= 2 4; A(R, A donde 4,, 


cualquier número tal que m,, < 1, , S M, ,, y límite de la misma cuando el diámetro 
dl ,q) = maxtd(p), d(a)) tiende a 0 a un número real L para el que se verifica que, 


para todo £>0 existe un 9>0 tal que, si d(p,q)<ó9, entonces IL-o(p, q) <E; 


decimos entonces que f' es integrable si y sólo si existe este límite, al que llamamos 
integral de fenR . 


Se demuestra que ambas definiciones son equivalentes, es decir, tanto las funciones 
integrables como los valores de sus integrales son los mismos. 


Integrabilidad de las funciones continuas 
Teorema. Toda función continua en un rectángulo es integrable en él. 
O, más general: 


Teorema. Toda función acotada, y continua excepto en los puntos de un número finito 
de curvas de clase C*, es integrable. 


Admitiremos estos teoremas sin demostración. 
Integrales iteradas 


Sea f una función definida en un rectángulo R como antes. Para cada x ela,b] queda 
definida en [c,d] una función f, de y, haciendo f, (y) = j (, y). Si esta función es 
integrable, tenemos la integral simple 


o =[ (y), 


que naturalmente es una función de x definida en [a,b]. Si ésta a su vez es integrable, 
tendremos la integral 


1 sy) dyjas. o simplemente DF riraar. 


a la que llamamos integral iterada de fen R. 


Cambiando el orden de las variables y con hipótesis análogas Antnirmos la otra 
posible integral iterada en R, 


ln es (x, y) dx dy . 


Ejemplo. Tomemos la función f la y) =x* +xy y el rectángulo R= [o, 1] Xx [- 2, 3] ; 
Tenemos 


13 1 xy” y=3 5% Sí? 5y? 1 35 
a 


Teorema de Fubini 
La integral doble se relaciona con las iteradas, de la siguiente forma: 


Teorema. Sea f una función integrable sobre un rectángulo R=|a,b]x[c,d], y 
supongamos que para todo x e[a,b] la función f,, definida como antes, es integrable 


en [c,d|; entonces, la función 


xo | r(xy)dy 


es integrable en [a, b] , y 


[40 9) dy dx = |, f(x, y) a 


Y lo mismo para la otra integral iterada, suponiendo ahora la integrabilidad de las 
funciones f, , definidas para cada y € [c,d] de manera análoga. 


Admitiremos este teorema sin demostración. No obstante, podemos tener una buena 
intuición geométrica del mismo. Si fes una función integrable positiva, su integral en R, 


que es el límite, para d(p,q) > 0, de la suma de los volúmenes de m-n prismas 


elementales de base En ] Xx | Y ¡19 Y ,] y altura 4, , (recordemos las sumas de áreas de 
rectángulos elementales en la definición de la integral simple), es, por definición, la 
medida del volumen cubierto por la gráfica de f sobre R, es decir, el volumen de 


2=((x,y.2)[a<x<b,c< y<d,052< f(x,y)) c R”. 
Por otro lado, para todo valor x, € [a,b] : 
[rw 


es el área A(x,) de la intersección de £2 con el plano x= x,. La integral iterada, 


[Acgar. 


no es entonces otra cosa que el límite de la suma de volúmenes elementales 


2 Al), .. Xu) 


A 


SOS 


NU: 


SS 


EOS 


ES 
OSI 


SS 


ISS 


SS 


== 


(cilindros de base dicha intersección, de área A(x,) , y de altura x, — x,_,: ver figura); es 


decir, es también claramente la medida del volumen de 42: esto es lo que afirma el 
llamado principio de Cavalieri (Italia, s. XVID), aunque es una idea que se encuentra ya 
en Arquímedes, s. III a. C. El mismo razonamiento vale para la otra integral iterada. 


Naturalmente, una de las cosas que afirma el teorema es que, si f es integrable, ambas 
integrales iteradas son iguales. Comprobémoslo en el ejemplo de más arriba. 


3 2 quel 
[ferojas- [52 To- [Do 


que es, efectivamente, el mismo valor obtenido antes. Esta permutabilidad del orden de 
las variables puede tener consecuencias importantes en la práctica del cálculo de : 
integrales, como veremos en el siguiente 


Ejemplo 


Sea f(x, y)=xsenxycos1rx, y R=[1, 2]x[0, 77]. 
2 px Tr 2 
[56 y)dA = f Jn x sen xy Cos Tx dy dx = f fx sen xy cos Tx dy dx, 
integrales iteradas iguales según el teorema. Pero la primera se calcula fácilmente: 
f fox sen xy cos Tx dy dx = f [- COS Xy COS rx), A = IM (- cos? TX +COS Tx) dx = 


1 


=f - A ox Js =>: 


mientras que en la segunda llegamos, tras aplicar identidades trigonométricas e integrar 
por partes, a 


2cos2y+eosy , (9 +7")ben2y+ emy) 


[ xsenxy cos zx ds =-217 a > a 


> 


y ¿cómo seguir? 
Propiedades de la integral doble en un rectángulo 


A partir del teorema de Fubini se extienden inmediatamente a la integral doble algunas 
propiedades básicas de la integral simple. 


1 - Aditividad respecto del dominio. Si R=R,UR,, siendo R, y R,dos rectángulos 


cuya intersección es de área nula (es decir, sólo tienen un lado común), y f es una 
función definida e integrable en R, entonces f'es integrable en R, y en R,, y se verifica 


que 
f aa = |, / da +), F da. 


2 - Linealidad. Si f y g son funciones integrables en un rectángulo R y k es un número 
real, entonces las funciónes f +g y kf lo son también, verificándose que 


[(r+g)aa=!, 7 aa+), da. 
| fda=k) 1 da. 
3 - Si fes integrable y positiva en R, entonces 
f fdaz0. 


Como consecuencia de 2 y 3: 


4 - Orden. Si f y g son integrables en R y f < g, entonces 
f, fda< |, gda 


(puesto que O < Jr -g)aa =/_ $ d4- |, £d4). 


Como corolario de lo anterior, y recordando que la integral de una constante en R es su 
producto por el área de R, resulta: 


5 - Teorema de la media. Si fes integrable en R y son « y [P respectivament su supremo 


y su ínfimo, entonces existe un número real u € [£, al, llamado media integral (o valor 
medio) de fen R, tal que 


| fd4=u0-ald—-c)=uAR); 
y sif'es continua, 1 será un valor de la función, es decir, 


[,£d4=y(Em an) 


para algún (£7)eR. 
Integral doble en un dominio no rectangular: definición y propiedades 


Hasta ahora hemos hablado de integrales en rectángulos, pero necesitamos ampliar el 
concepto a dominios más generales. Y lo haremos limitándonos en lo que sigue al caso 
de funciones continuas, quizá con la salvedad indicada en el segundo teorema de la 
página 3. 


Sea D un dominio del plano limitado por una curva C de clase C* e incluido en un 
rectángulo R, y sea f una función continua en D. Definimos en R otra función, f,, de 


la siguiente forma: 


O enotro caso. 


pda ¿E si (x,y)eD, 


La función f, es obviamente continua excepto en una curva de clase C*, luego es 
integrable en R. Entonces, por definición: 


[da = [fo da. 


Las propiedades de la integral sobre un rectángulo que veíamos más arriba: aditividad 
respecto del dominio, linealidad, positividad, orden, así como el teorema de la media, 
se extienden de manera natural a la integral doble sobre estos dominios más generales. 
Y, si f es positiva, su integral admitirá la misma interpretación geométrica como 
volumen del sólido 


Q= (x, y, 2) (x, y)eD,0<z< f(x, y»). 
Y tomando f(x, y)=1, será también 
[ da =A(D), 


es decir, la integral doble será la medida del área del dominio D. 
Dominios elementales del plano 
Llamaremos: 


región elemental de tipo I a una parte del plano que 
pueda definirse así: 


D=4(x,y)]la<x<b, o) <y<p,0), 


; donde , y f, son funciones de clase C ' en [a,b], la 
a primera mayorada por la segunda; 


región elemental de tipo II a una que pueda definirse como 


D=4(x,y)|e< y<a, yw (y)sx<w,(»)). 
con las mismas condiciones para Y, y Y,; 


región elemental de tipo III a una región que pueda definirse lo mismo como de tipo 1 
que como de tipo II; y finalmente 


dominio elemental a una parte del plano que sea unión finita de regiones elementales de 
cualesquiera de los tres tipos disjuntas dos a dos, o bien con intersecciones de área nula 
(formadas sólo por puntos frontera). 


Integrales doble e iteradas en dominios elementales 


Sea f una función continua en una región elemental D de tipo L y c y d cotas, 
respectivamente inferior y superior, de y, y de Y,. Si es R el rectángulo [a,b] x[c,d] y 
definimos en él la función f,, como antes, será: 


J,rda=f, 1 da= 5 Lr,dyax= (1% 5, dy + 


p,(x) 


roto lt) das, 


ya que las integrales primera y tercera dentro del 
paréntesis son nulas, por serlo f,, en esos intervalos, y en 


la segunda f,, y fson idénticas. 


a b 


Análogamente, si D es de tipo Il, se verificará que 
vw, (») 
f rds= legis fdxdy, 


con demostración idéntica. 


Si D es de tipo III, la integral de fpodrá expresarse como integral iterada de dos formas 
distintas: 


e d qwa() 
d, Saa al dy ds =| de Fdxdy. 


V q) + 


Y finalmente, si se trata de un dominio elemental D = U D, , con las regiones D, dos a 
i=] 


dos disjuntas o con intersecciones de área nula, 


|, sas = Y), sad ] 


Ejemplo 1 


Sea f (x, y) = xy y D la región del plano encerrada por las curvas 


y=6-5, y=2WHx, x=0. 


La intersección de las dos primeras tiene de abscisa x=4, y 


Xx 
entre este valor y x=0 es da >2/x, luego 


p=((y)losxss 3 y<6-3), 


con lo que 


Boeoroooroooroo=oo.- 


Ejemplo 2 
Sean ahora f (x, y) =0” y D el triángulo determinado por las rectas 
x=0, y=2 e y=2x, 
es decir, 
D=Á(x,y)0<x<1, 2x<y<2). 
La integral será 
|raa=[[ “aa: 
D 02; dió 


que no podemos resolver por no existir una primitiva elemental de 


Pero D es una región elemental de tipo TI, que podemos escribir 
también así: 


0sys2.05x5%): 


D= (5) 


con lo que 


y=l Y 
A E 
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lo que ilustra la importancia práctica que puede tener el orden de las variables en la 
integración iterada. 


Geometría de las aplicaciones de R? en R? 
Para abordar el importante tema del cambio de variables en una integral doble tenemos 
que hablar de aplicaciones del plano real en sí mismo, y de cómo éstas afectan al área de 


un dominio. 


Sabemos que, si 9 es una aplicación lineal e inversible (isomorfismo) de R? en R?, 


x = qu + bv la É a bYu 
y=cu+dv " e y) lc div 


donde 


es una matriz regular, entonces 0 transforma rectas en rectas y conserva el paralelismo, 
y además el área de un paralelogramo se ve multiplicada, en su transformado, por el 
valor absoluto del determinante de 4. 


Sea por ejemplo la aplicación 


u 3v 
x==+— 
2 L 
y =-u+2v 
Tenemos: 
Y 3 
|det A = Al 4 ma 
ii 4 2 
2 
Sea D' el cuadrado [o, 1]x[0, 1], de área 1; como 
1.5 


: Ñ 50, 0)= (0,0), DÍ, (3-1) 
(5.2) D(1, 1)= (2, 1), 


g(0, 1) 


-0.5 . , re 
D” se transforma en el paralelogramo D de vértices 


dd (0, 0), (> - .) (2 2) y (2,1), cuya área es 2, como 


se comprueba sin dificultad. 
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Si la aplicación D no es lineal, pero sí de clase C* e inversible en un dominio D, lo 
anterior no es inmediatamente aplicable. Pero sabemos que existe entonces en cada 
punto, localmente, una buena aproximación lineal, a la que llamamos diferencial de D, 
cuya matriz es la de derivadas de DP o matriz jacobiana, y cuyo determinante, el 
jacobiano, no se anula por ser Y biyectiva. Naturalmente, este determinante ahora no es 
constante, depende del punto. Pero localmente nos permite relacionar áreas de partes 
muy pequeñas del plano con sus transformadas por O. 


Veamos un ejemplo. Sea la aplicación DY definida así: 


x=u +uv 
y=vw -u 
cuyo jacobiano es, en cada punto, 
A 2u+v uy 
y 4229) - = duv+2v?+u. 
(u, v) =1 2 


Consideremos como antes el cuadrado 
D' =[0, 1]x[o, 1]; 


su área es A(D”) =1, y el valor del jacobiano en un punto interior, por ejemplo el 


centro, 


con lo que 


Hallemos ahora el transformado D = o(D'). Los vértices se transforman en 


D(0,0) =(0,0), D(L0) = (1,-1), 
D(1,1) =(2,0), 09(0,1)=(0,1); 


pero ahora las rectas ya no se transforman en rectas, sino que el lado v=0, por 
ejemplo, se transforma en 

n= 8 

Y == 


es decir, en la parábola 
y en Lx > 


y análogamente los otros tres lados, según el cuadro 


u=lb y=x?*-2 


3-y1+4x 
2 


v=lhb y= 


u=0hbx=0 
con lo cual D es como indica la figura; su área es 


A(D) = [A [Vx as [ (1 -2x) dx = 


2 3 1 


ax (ra ] [28] e 2 
2 3 gl 
0 0 
23 28-42 41)=8 -21666...; 
Bus. 3 6 


la aproximación es evidente, y sería mayor para dominios D" más pequeños. 
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Esto sugiere lo siguiente: puesto que la integral es el límite de la suma de áreas muy 


pequeñas del plano, multiplicadas por un valor de la función en un punto de ellas, 
Lr(xy)axay = lim 4 (x.y) 45,4, , 
¿no habría quizá que sustituir, al cambiar de variables mediante la aplicación 
9:D">D, 
el factor Ax, Ay, por o (u, y, ) Au, Av, ? Es decir, llamando 


F"(u,v)= £(x(u,v), y(u,v)), 


¿no será 


[£(x,y)axay = lim $ "(uy vi) (our) 40, 4, = [..£"(u,v)Vo]dudv? 


Tomemos un ejemplo sencillo. Sea D el círculo x? + y? < a?; sabemos que su área es 
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[aa a. 
Hagamos en esta integral un cambio a coordenadas polares, 
(x, y) = Dlr, 0) = (7 cos 0, r sen 0) a 
cuyo jacobiano es 


coso —rsen0 


sen0  rcos0 


Evidentemente 
D= o(D'), siendo D'= ((r,0)| 0O<r<a,0<0< 27) 


(0 es de clase C* e inversible en el interior del rectángulo D*, donde r no se anula), 
con lo que tendremos 


[rarao=[| [Fardo= 27 | rdr= | a 
0 


es decir, en este caso sí se verifica la conjetura. 


Pues bien, ello no es casual, como afirma el teorema siguiente (del que todo lo dicho no 
es, desde luego, una demostración, pero sí una aproximación que lo hace intuitivamente 
más natural). 


La fórmula del cambio de variables 


Teorema. Sea D un dominio elemental del plano, y sean DP: D" > D una biyección de 
clase C* y funa función continua en D (pudiendo no verificarse estas propiedades de 
9 y de fen un número finito de curvas C*). Entonces 


[120-Vol=[£; 


o dicho de otro modo, 


EA dy 
, dudv= | f(x, y) dx dy. 


(u,v) 


J 


f. f(x(u, v), y(u,v)): 


Ejemplo 1 


Vamos a calcular el área del dominio D definido por las curvas 


; y=x”, iaa”, we=l, xp-8 
en el primer cuadrante. 


Tomaremos como nuevas variables 


00.51 1.5 2 2.5 3 


con lo cual el dominio se convierte en el rectángulo 


(1, y) 


(u,v) 


y el jacobiano J será el inverso de 


J (a, y) e E 3 = Ay 
An Xx e AE Ta 
(x, y) $. % Y? 
es decir, en valor absoluto y puesto que uy es positiva, 
(x, y) 0 O 
AAA, 
(u, v) 4y| 4u 


con lo que resulta finalmente 


1 4 
7 


l '  7log7 
en . EL 234053. 


dv du = que 


A(D)=Jif, . 


1 
4u 


Comprobémoslo calculando el área sin cambiar de variables. Las intersecciones son: 


y=x* = y=x z 
= y== = A 
1,x=l; 2 =07; 3.2108; n= 4/56; 
LE y=-= LT y=-= 

Xx Xx 
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-7-1_187,68-7 ¿log56-log8 _56-8 _ 


4 4 7 4 4 28 


42+6-48 _ 7log7 
28 4" 


Ejemplo 2 


Sea D de nuevo el círculo x? + y? <a?, y f(x, y) = x?y?. Calculemos la integral de f 
en D mediante un cambio, de nuevo, a coordenadas polares. Tenemos: 


f"(r,0) = f(o(r,0)) =r*cos* Osen* 0, 


con lo que 


[ss yaa = [ fr cos* O sen? OdO dr =[ a e ar = 


F. 0=27 6 ]1 6 
=[[ E rro E A ] 
8 S 3% les , 


Aplicaciones de la integral doble 

La integral doble admite aplicaciones geométricas y físicas importantes, como son la 
definición, cálculo y propiedades de algunos conceptos y magnitudes básicos, tales 
como volúmenes, masas, centros de gravedad y momentos de inercia. En lo que resta 
del tema se hablará de todo ello. 


Volúmenes 


Ya hemos visto que la integral de una función positiva f en un dominio D mide un 
volumen: 


[raa=v(0), 
donde 
52 = ((6.9,2)| (x, y) ED,0<z< £x,y)). 


Ejemplo 1 


Sabemos que, al cortarse ortogonalmente, dos bóvedas de medio cañón de la misma luz 
engendran dos bóvedas de esquina, una formada por la hoja superior, bóveda de arista, 
y otra por la inferior, bóveda de rincón de claustro (ver figura). 
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Calculemos el volumen cubierto por una de ellas, 
la de arista por ejemplo, suponiendo que la luz de 
las dos de medio cañón sea 2/ metros. Tomando 
ejes de coordenadas de la manera obvia, las 
ecuaciones de los dos cilindros serán, 
respectivamente, 


cerrar, Pita, 


y un octavo de la bóveda de arista quedará definido, por ejemplo (ver plantas), por el 
cilindro z = 1? — x? sobre el triángulo D encerrado, en el plano OXY, por las rectas 


To 
=== es decir, el volumen total cubierto por la 
A bóveda será 


yA y=x,y=l, x=0, 


V=8/ 2dA, 


x=lsent, dx=1lcost dt, 


8 7 
Vi? —-x? =lcost; x=051=0, =let=>; 
a V =8 [*1"(1- sent) cos? tdt = 


“(1+c082t 
8p PE co sent) = 


n 


S. 
= (27-5) *m?*=3616521* m*, 


0 


sen2t  cos?t 
+ 


=81| 24 
2 4 3 
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lo que supone el 90,41% del correspondiente a la cubierta plana de igual altura, 4/* m?, 


De forma semejante, cambiando, bien la superficie, bien el dominio de integración, 
puede calcularse el volumen cubierto por la bóveda de rincón de claustro. 


Ejemplo 2 


Calculemos ahora el volumen del menor de los 
segmentos esféricos definidos por la esfera 


+ y? +2? =2]? 
y el plano 


z=k, 


La proyección de la intersección de ambos sobre el 
plano OXY será la circunferencia 


y el volumen buscado 
V = [zaa , 


donde 


z= 4/21? —- x? - y? =1, D=4(x, y)|x? + y? <P). 


Cambiando a coordenadas polares, resulta 


¿=V42P -r? -1, D'=((r,0)0<r<1,0<0<2n), 


V= [ (V2F=7 -1)raodr= 2 [lar y? +11) dr = 


19 


A partir de este volumen, y recordando que el de una 


4 
esfera de radio R es ZAR, es inmediato calcular el 


cubierto por una bóveda baída, definida por una 
semiesfera sobre un cuadrado de 2/ metros de lado, es 
decir, apoyada en cuatro arcos de luz 2/ m. (¿cómo se 
haría? se sugiere comenzar considerando la esfera 
completa, tal como representa la figura; para 
comprobación del que aborde el cálculo, diremos que el 
volumen pedido equivale al 80,40% del encerrado por 
una cubierta plana de iguales planta y altura). Y el 
siguiente paso sería calcular el de una bóveda 
semiesférica sobre pechinas esféricas. 


Supongamos ahora que queremos calcular el volumen de 
un sólido, 2, definido por dos superficies dadas 
explícitamente por dos funciones continuas f y g, 


Z = f(x, y) 
Za = g(x,y), 


que se cortan según una curva regular C, cuya proyección sobre el plano OXY sea una 
curva plana simple C, que encerrará un dominio D. Es claro que el volumen definido 


por dichas superficies será 
v(2)= [|s- glaa 


(de forma análoga a lo que ocurría, en la integral simple, con las áreas de figuras planas 
encerradas por dos curvas). 


Ejemplo 1 
Supongamos que queremos calcular el volumen encerrado por el paraboloide 
2=x* +3y? 
y el elipsoide 
2x7 +6y+2?*=3 
en el subespacio z > 0 (ver figura). 
La intersección de ambas cuádricas estará en el plano dado por la ecuación 


22+2"=3, obien (2-1)(z+3)=0, 
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con la condición z 20, es decir, z=1, que llevado 
a la ecuación de cualquiera de las cuádricas nos da 


x4+3y =1, 


elipse de semiejes 1 y 5 que define, en el plano 


OXY, el dominio D de integración. Y como en dicho 
dominio el elipsoide está por encima del 
paraboloide, será finalmente: 


y = [ (3-27 =6y* -x* -3y*)aA. 


Para resolver esta integral haremos el siguiente cambio de variables: 


x=rcos0 
1 
ll 


(parecido al cambio a coordenadas polares, pero que transforma rectángulos en elipses 
y no en círculos como aquél; se habla por ello, a veces, de coordenadas elípticas), cuyo 
jacobiano es 


o 
ES 
< 

_— 


r 


y ; 
(r, 0) 3 


con él, el dominio y las dos superficies se convierten respectivamente en 


D'=((r, O l0<r<10<0<2%), z2=r* y 2=v3-27, 


resultando 


r=0 


Si la intersección de las superficies no fuese una curva simple, podría ocurrir que en 
unas partes del dominio una superficie estuviese por encima de la otra y en otras por 
debajo, con lo cual habría que calcular varias integrales cambiando el signo, etc. 
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Pero el cálculo de volúmenes puede abordarse también mediante la integral triple, por lo 
que aquí no diremos más de ello. 


Masas 


Si en un dominio D del plano tenemos una función positiva, podemos interpretarla 
siempre como una densidad, medida naturalmente en unidades de masa por unidad de 
superficie. La integral medirá entonces la masa del dominio: 


M(D) = [mx y) aa, 


donde m es la función densidad a la que nos referíamos. 


Cualquiera de los casos de integrales de funciones positivas que hemos visto podría 
servir de ejemplo (en el siguiente apartado veremos uno específico). 


Centros de gravedad 


Dado un sistema de n masas m,,m,,...,m, , situadas sobre una recta en los puntos de 

abscisas respectivas X,,xX,,...,X,, sabemos que se llama centro de gravedad, 0 

baricentro, del sistema al punto respecto del cual el momento (estático o de primer 
n 

orden) del sistema es nulo, es decir, el punto de abscisa x tal que Nx, x) m, =0, 
k=1 

de donde resulta 


n 
donde M = Nm, es la masa total. 
k=1 


Esto se generaliza al caso de dos dimensiones de la forma obvia, y a un continuo de 
masas, es decir, a un dominio D del plano con una función de densidad m:D>R,, 
mediante la sustitución de sumas por integrales. En resumen, las coordenadas del centro 
de gravedad de D serán 


7 trás jaa _ [ema y)an _  ymlsy)aa 
o | m(x, y) da Ñ M y PA ; 
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y estaríamos ante un concepto puramente geométrico. 
Ejemplo 1 


Calculemos el c.d.g. de un semicírculo de radio a. Tomando el centro del círculo como 
origen de coordenadas y el diámetro frontera como eje OX, será x=0 por simetría, y 
con un cambio a polares, 


meta sen0d0 dr af” 2 dr 4a 


Ta Ta =3p 


Es decir, el centro de gravedad es el punto situado sobre el eje de simetría a una 


4 
distancia pS = 0,4244 a del centro. 


Ejemplo 2 


Calculemos ahora el c.d.g. del mismo semicírculo pero con densidad igual, respecto de 
la misma referencia, a k (a * y) , donde k es una constante dada. 


Al no variar la densidad con la coordenada x la simetría se conserva, con lo que seguirá 
siendo Xx = 0. Calculemos, pues, primero la masa y luego la coordenada y : 


O=x 
M= [| [r(a=rsen0)rd0dr = [rlao+ rcosó| dr= 
0=0 


2 ae 
- [tr(ma—arjar=k ES = | E a 
0 


My= [ [r(a=rsen0)r? sen0d0 dr = [ [ tr? (a seno- 


1-cos20 0 1-8 
a [rr[ =acoso- 70% | dr= 
2 2 4 Jozo 


3 471 e 
= [tr 2 dre 2ar?  7r ¿16 3 
31” a 


16-37 , 6 16-37 
a A AN ETA 
24 k(31-4)a 
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Momentos de inercia 


Es sabido que el momento de inercia (o momento de segundo orden) de un sistema de n 
puntos, de masas m,, m,,...,m, , respecto de un eje e (magnitud física que desempeña 


n 
ante un par de fuerzas el mismo papel que la masa ante una fuerza), es /, = y d'm, , 
i=1 
donde d, es la distancia de la masa m; al eje e. Esto se extiende a un dominio D del 
plano, de densidad m: D > R, , mediante la integración. 


Si tomamos el eje e como uno de los ejes coordenados y es míx, y) la densidad en cada 


punto, el momento será 


I, = [mx y) y dA, 


o bien 
de = [mr(x,y)x" dA, 


según que e sea el eje de abscisas o el de ordenadas. 


También se habla de momento de inercia respecto de un punto, o momento polar, que 
equivaldría al momento respecto de un eje perpendicular al plano. Si el punto es el 
origen, evidentemente 


l)=1,+1,. 


Ejemplo 1 


Calculemos el momento de inercia de un rectángulo, de lados a y b y densidad constante 
m, respecto de su base media paralela al lado de longitud a: 


b 
. 312 by DM 
1,= [my ardid El > A 
2 me 
2 


Ejemplo 2 


Sea ahora el semicírculo del ejemplo 2 de la página 21. Calculemos el momento de 
inercia respecto del eje de simetría: 


1+c0s20 
1, = [ [xa -rsen0)r? cos* OrdO dr = k[ a 
A 39 O=x 
— r?* senOcos* 0)d0 dr = kl arq +y* ex] dr = 
0=0 
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3 4 
rar” 2 Tm 2 
ef A (4-23) a. 
2 3 8 15 

El siguiente teorema enuncia una sencilla propiedad relativa a los momentos de inercia 
respecto de dos ejes paralelos. 
Teorema. El momento de inercia respecto de un eje cualquiera es igual al momento 
respecto del eje paralelo que pasa por el centro de gravedad más el producto de la masa 
por el cuadrado de la distancia entre ambos ejes. 


Demostración 


Basta ver que, si el centro de gravedad es el origen de coordenadas y el eje dado la recta 
y =d,es 


(y-d) = y -2dy+d”, 
es decir, 
[, = [6-0 m(x, y) dA = La -2d | y m(x, y) da + Md? = 1, +Mad?, 


ya que la integral del segundo miembro es el momento estático respecto de un eje que 
pasa por el centro de gravedad, y por lo tanto es nula. M 


Por ejemplo, ¿cuál sería el momento de inercia del rectángulo del ejemplo 1 de la página 
23 respecto de su base? Según el teorema, 


ll - Integral triple 


Integral triple en un rectángulo 
Todo lo que se dijo para introducir el concepto de integral doble sirve, con las pequeñas 
adaptaciones necesarias para pasar de dos a tres variables, como introducción de la 


integral triple. Ahora consideraremos una función acotada 


fiR>oR, 
donde R es un rectángulo (paralelepípedo rectángulo) en R?, es decir, 
R=|a,, b,]x[a,, b, ]x[a,, 5,], 


y sendas particiones p, q y r de los tres intervalos a, Bl, que dividirán R en m-n- p 


pe 7 


rectángulos elementales R,,,, en cada uno de los cuales f tendrá un ínfimo y un 


supremo, M, ,, y M,,, respectivamente, con los que definiremos las sumas s(p, q,r) = 
= 2, Min v(R,,,) y s(p,q,r)= a v(R, ,,), donde VR.) es el volumen del 
Dis DI 


rectángulo elemental R,,,, etc., resultando la integral como punto frontera entre el 


conjunto de las primeras y el de las segundas, si es que existe, que será la condición para 
que f sea integrable. 


También la definición como límite es análoga ahora, a partir del concepto de suma de 
Riemann, p(p,q,r)= A v(R, ,.), donde Ju, ,, es cualquier número real tal que 
ij k 


Mx SM, ¡x S M, yy, y de su límite cuando d(p,q,r)= maxld(p), d(g), d(r)| tiende a 


0. Ambas definiciones son equivalentes. 


Las notaciones para la integral triple son 
1,.£ | fas dy dz, 1 y, 2) dx dy dz, SJ, se y, 2) dx dy dz, 


o bien escribiendo dV en lugar de dx dy dz. 


Integrabilidad de las funciones continuas 


Teorema. Toda función definida en un rectángulo R, acotada, y continua excepto en los 
puntos de un número finito de superficies de clase C*, es integrable en R. 
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Integrales iteradas 


Serán ahora de la forma 
bi pb, Pb, 
L lA IN F(x,y,z)dzdydk, 


supuesta la integrabilidad de todas las funciones que se van integrando sucesivamente 

(las f, , definidas por f, ,(z) = flx,y,z), ze [a,, 5, ], y las f, definidas por f,(y)= 
b 

= j 3 de (x, y,z)dz, y€ [a,,b, |), condición que se cumplirá si f es continua, pues ello 


implicará la continuidad de todas esas funciones. 


Naturalmente existirán hasta seis integrales iteradas distintas, pues seis son las 
permutaciones posibles de las tres variables. 


Teorema de Fubini 


Afirma, igual que en la integral doble y con análogas condiciones, la igualdad de la 
integral triple y de las seis integrales iteradas. 
Ejemplo 


2 


Sea f la función f (x, y, z) = = +z,y RCER? el rectángulo [0,2]x [1,4]x [- 1,1]. Será 


S, 46 y, z)dV = ¡N E +2) dy dz dx = f E [7 l0g y +2, dz dx = 


2 


=/f Qlog2x” +30) = [ligar dx =|4log2x” dx = mt 
0 9-1 0 2 Local 0 


3 


Integral triple en un dominio no rectangular: definición y propiedades 


La definición la haremos como en la integral doble, extendiendo la función f mediante 
valores nulos a un rectángulo R que contenga el dominio £2, y haciendo entonces 


[_rav=5£, av. 


Con esta definición la integral triple tiene las mismas propiedades que veíamos para la 
doble, es decir, aditividad respecto del dominio, linealidad, positividad, orden, y se 
verificará también el teorema de la media, es decir, existirá un real 41, acotado por el 
ínfimo y el supremo de fen £2, tal que 


[17 =uv(o), 
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siendo 4 = f (E 7, c) para algún (E 7, 5) e 42 si fes continua. 


A partir de ahora consideraremos, en general, funciones continuas (con la posible 
salvedad ya hecha en el teorema sobre la integrabilidad de éstas), y hablaremos también 
de regiones y dominios elementales. 


Una región elemental de tipo Tes toda parte del espacio que pueda definirse así: 


albacea), 


siendo D un dominio elemental del plano OXY, 
y YW, y Y, funciones de clase C* en D. Y 
regiones elementales de tipo II o tipo III serán 
las que puedan definirse cambiando el orden x, 
y, z de las variables a z, x, y O y, Zz, x 
respectivamente, llamándose región elemental 
de tipo IV la que sea a la vez de los tres tipos 
anteriores. Y finalmente, un dominio elemental 
será la unión de un número finito de regiones 
elementales disjuntas dos a dos, o bien con 
intersecciones de volumen nulo (sólo puntos 
frontera). 


Si 42 es una región elemental de tipo ly D= dy) as<x<b, o (x)<y<0, 07 será 


l f (x, y, 2)dV = ÍN 17 lentes f (x, y,2)dz dy dx, 


ac) y 


y lo mismo en los otros casos, permutando las variables. Y si se trata de una región 


n 
elemental, £2=1742, con las condiciones sabidas para las regiones £2,, será 


i=1 


|,147 Y 5,547. 


Ejemplo 1 


Sea (2 = X(x, y,2)| 0O<Sxs10<y<l-x,0<zS1l-x -y) (¿qué figura geométrica es, 


y qué tipo de región elemental?) y f (x, y,z) = xyz . La integral será: 


z=l-x- 
ron] a 


z=0 
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2 4 


2x[l-x)y? aye pa" 11 1 > 4 
PE z Y — e dx = | A x(1-x)" dx = 


y=0 


1 | Ea El A 
“Ao ro. 6 | 720 


Ejemplo 2 


Calculemos el valor de 


z 


Á4= [LJ dz dy dx 


mediante una integral iterada. Por ejemplo, 


4=[[72 e” dl a 


pa 


x = [5 De” dydx=.. 


no podemos continuar; pero el dominio es una región de tipo IV, luego podemos elegir 
como queramos el orden de las variables: por ejemplo, 


2=4(x,y,2)]0<x<1x<y<1y<zS1)= 
=((x, y,2)J0<7<1,0< y<z,0<x< y); 


con lo cual 


4 [ff ada [aya [230 ¡E 2 => [ee -2]dz= 


== [e -2*] =3(0-1-1)= e 


lo que ilustra aquí de nuevo la importancia que puede tener el orden de las variables al 
integrar iteradamente. 
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Cambio de variables 


Todo lo relativo a la geometría de las transformaciones en R? es aplicable punto por 
punto en R?, cambiando ahora áreas por volúmenes, con lo que, recordando las 
conclusiones a las que allí llegábamos podemos enunciar: 


Teorema. Sean 42 un dominio elemental del espacio real, f una función continua (con 
las salvedades conocidas) en 42 y D:22" >.2 una biyección de clase C'. Entonces 


[sf = [ £-0-Vol, 
o bien 
IE y,2)dV = l. qu) E22) du dv dw. 
Ejemplo 1 
Calculemos la integral de la función continua f:22>R, donde f(x, y,z) = een y 


Q= ¿(, y,2)| 1<x* + y <4)1<z< 2) , mediante el cambio de variables 


x=rcos0 
y =rsenU0 
LZ=Z, 


cuyo jacobiano es 


(que sólo se anula en una recta, el eje. OZ, por lo que no hay problemas con la 
biyectividad). Será: 


f"(r,0,2)=1(0(r,0,2))=%, 2" =Í(r,0,2)]1<r<2,0<0<2x,1<2<2), 


[, lu. .Jav = [, 5(7,0,2)rardodz= | | ["Hzrdodzar = 


=27 | Parar=an [E ara [2ar=3m0g2. 
z=l 
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Naturalmente, este cambio de variables ha sido a las que conocemos como coordenadas 
cilíndricas, que consisten simplemente en tomar coordenadas polares en el plano OXY 
dejando la z invariable (la ecuación r = k es la de un cilindro de revolución, de radio k 
y eje el OZ). 


Recordemos igualmente las coordenadas esféricas, definidas por las relaciones 


x=pcosOsenq 
y=psenOsenp 
Z=pcosQ 


(p es ahora la distancia del punto al origen, y p =k es la ecuación de una esfera de 
radio k y centro el origen), cuyo jacobiano será 


(x, y,2) cosOsenp —psenOsenp pcosOcosq 


Ll 
(p,0,p) 


=|sen0senp  pcosOsenp  psenOcosqp = 
cos p 0 = psenq 


=p? (cos p|- senpcos p) — sen p sen? p) =-—p? senq, 
que se anula también sólo en el eje OZ. 


Ejemplo 2 


Calculemos la integral de f(x, y, 2) =wyx? +y?+2? en la esfera de centro el origen y 
radio a, £2: 


2 =((p,0,9)/0<p<a,0<0<2x,0<p<r), F (p0,0,p) =p, 


a pr p2x a =7T 41" 
14 dV = f f f p" senpdOdpdp = 27 pr cos pl, dp= an E =a!, 


Aplicaciones de la integral triple 

Lo mismo que en el caso de la integral doble, las aplicaciones de la triple serán 
geométricas, en el cálculo de volúmenes, y físicas en el de masas, centros de gravedad y 
momentos de inercia, ahora de sólidos en el espacio de tres dimensiones. 


Volúmenes 


Naturalmente, dado un recinto 42 tendremos simplemente: 


v(2) = |, dv. 
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Ejemplo 1 


Calculemos el volumen de un elipsoide de semiejes a, b y c. Para ello haremos el 
cambio obvio 


x=apcosOsenq 


y =bpsen0senq 


Z=CPCosQ, 
cuyo jacobiano es 
lx, y, z) 3 
JFK =-—abcp” senq 
(p,0,9) 


y con el que el elipsoide referido a sus ejes se convierte en 
p=1, 


siendo por lo tanto su volumen 
lp px 1 E 4 
V(<2)= NN ' abcp? senp d0 dp dp =27 abc | p* [- cos ol do = zFabe 


Ejemplo 2 


Sea ahora el sólido 2 definido, en coordenadas 
esféricas, por la superficie 


0=0p 
y la esfera 


p=RkR 


(ver figura). En virtud de las evidentes simetrías y 
para conseguir que el signo de sen, y por lo tanto 
el del jacobiano, se conserve positivo (lo que 
facilitará el cálculo del volumen), consideraremos 
sólo un cuarto del sólido: 


h 


> 


SS 


SS 


SS 
O 
SS 


SNS 
XS 
SS 


TT 


,0<gp<Q0»; 
2 p 


Q, =3(p,0,p]0<p<R,0<0< 


con ello, el volumen será 
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v(32) = af f (o sen p dp d0 dp = 4[ El- p cosg|.,d0dp = 


z 2(77-2)R* 
=4[ f01-coso)aoap= [o [2-1) qo= EE. 


Ejemplo 3 


Consideremos una bóveda engendrada por un 
arco de medio punto que gira alrededor de uno 
de sus apoyos según un ángulo 0, a la vez que 
su luz crece proporcionalmente a 0, como 
muestra la figura. 


Si la luz de cada arco es k0, donde k es una 
constante dada, la ecuación de la superficie, 
en coordenadas esféricas, será 


p=k0senq, 


y el espacio cubierto, suponiendo que O varía entre 7 y 
377, 


2-1(05,) 


05 pSk0snp.1s0s3050<h): 


con todo ello, el volumen será 


v(2)= f [770 senpdododp= É [CELS sen pao dy = 


302 4 4 By . 
Es sen pag= DE e e] prnl 


3 4 3 2 
Sk77* dl L1.coste) 
, 1-2c092p++£0% | qq = 
- Pp cos 2p . p 
dis z | 
_Sz sep cd SE 
312 e 
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Caso de un sólido de revolución 


Sea f una función de una variable, definida y continua 
en un intervalo [a,b]. Supongamos que, en R?, 
hacemos girar la curva y= f (x) alrededor del eje OX: 


se engendrará un cierto sólido de revolución, 2. 
Vamos a calcular su volumen pasando a las variables 


a 
y=rcos0 
z=rsen0 


es decir, a unas coordenadas cilíndricas cuyo eje 
r=0 es el OX, cambio cuyo jacobiano es 
evidentemente r (en la figura, para mayor 


claridad, O varía sólo de O a Y La curva 


generatriz de la superficie que envuelve £2 será 
ahora r = IA (x) en el plano 0 =0, y el sólido 2 


será el transformado de 
2" =(x,r,0)Ja<x<b,0<rs|f(x),0<0<27]. 


Con lo que el volumen será 


2 


y=([0 ["raoaras=2% [ [Urarac=2n [| Oc [ ds 
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2 r=0 
expresión conocida que permite el cálculo del volumen mediante una integral simple. 
Ejemplo 


Calculemos el volumen del sólido engendrado por un arco de sinusoide entre dos ceros 
consecutivos. Será 


= dd 2 
a a e sa. 
2 2 4 


Masas y centros de gravedad 


Al pasar de dominios bidimensionales a tridimensionales las integrales serán triples y no 
dobles, pero los planteamientos, definiciones y cálculos serán en todo análogos a los 
vistos en el primer caso. Veamos algunos ejemplos. 
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Ejemplo 1 
Consideremos el sólido 42 limitado por las superficies 


z, =9x+14y-20 
z, =x* +4 y? +5x-—18y +32 


y cuya densidad en cada punto es 


mx, y,2) = Qi 


S 


calculemos su densidad media, es decir, m = y" 


La proyección sobre el plano OXY de la intersección del plano y el paraboloide elíptico 
será, haciendo z, = z,, la elipse 


x' -4x+dy?-32y+52=0, obien (x-2)? +(2y-8)' -16=0, 


o finalmente 


Esto sugiere el cambio 


x=2+4rc0s0 
y =4+2rsen0 
Z=Z, 


cuyo jacobiano es 8r. 


Las superficies son, después de operar, 


z, =54+36rc050+28r sen0 
z, = 38+ 36rc050+28rsen0 + 16r?, 


y el sólido 
2 =4(r,0,)J0<r<1,0<0<2x,z,$2<2,). 


Con todo ello podemos ya calcular el volumen y la masa (abreviaremos algunas 
simplificaciones): 
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V= [dv =[ [7 f srazaoar=[ (8r(z, -2,)d04r = 
=8[ [" r(16+16r?)404r =2567 [ (1 +r")ar = (2+2) 2567 =19271 <= 603,1858 ; 


2 
2 abr = 


id / [ (A ado dr = d [ [lar+2r sen 0) z 
- [(lar+2r sen0)64(1—r?)(19+2r” +9rcos0+7rsen0) d0 dr = 
pr IN l- 1?) )(n6r+8r +(66r? +4r *)sen0+36r? cosO + 


+18r* senOcos0 + Tr (1—- cos20)) d0 dr = E ((1- r?)(76r + 151*)27 dr = 


_5127 52% (33 61 5). 1057 
y 


2357 4 2 


76r - 61r? -15r*) dr = e 


= 1302,8813; 


y finalmente 


_ 103687 1 5354 
25 197 25 


Ejemplo 2 


Calculemos la masa y el c.d.g. de una esfera de radio a y cuya densidad es, respecto de 
un sistema de referencia con centro en el de la esfera, m(x, y, 2) =2a-Zz. 


M= [ [ Fla - pcosp)o? senpdO0dp dp = 


4 


2, 79=2 
8 
o 2 [|- 2ap* cosp- p? ne. do = 27 | 4ap* dp = E 
p=0 


Por simetría, 


x=y=0; 


MZ= [ [ [Focosp(2a = pcosp)o* senpdO0dp dp = 27 | [ap cos Q sen — 


4% a 3 2 4 cos? y á El 
=p" cos psenp)dpdp =2x | ap” sen” p+p 3 do = 
e p=0 
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Ejemplo 3 


semiespacio x > 0 por las superficies 


! 
l) *+y3-2=0 
iN cat 
INN E 
NN 2x-y=0 
/ A ¿=0, 


Para calcular su baricentro haremos el cambio de variables 


X= 
ll 


Nn = 
HN 
SS 3 Efe 


cuyo jacobiano vamos a calcular: como 


LA 
LOA E E 
(x, y,2) 6 0 Xx 
invirtiendo, 
(ey2_[x]_ 1 
(u,v,z)| |2y|  2|u| 


Por otro lado, 


Sea ahora el sólido 42, de densidad constante, limitado en el 
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con lo que 


1 
Isus2.0svsoses (uso 


o luna 

Y [E zara [Gr rea; (Ja 
lib) 

vz [LEO fe zgtcin aa | (lis) varas 

loz ra 

2- 


yy fp >, ed lus Ddvdu= 


2 


30 so >? 


1 


1 2 
“520 


Lp gua 2)» Xamda=l [A : += 


Vz= AZ aca (+) lr (14242) au 


q A Ll2+21 A. 1,2. AN 
aa Sl a a gp 96 


_ 15+16log 2 
E” TS a 
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El centro de gravedad es el punto 


8(11/2-8) 8(/2-4) 15+1610g2 
a A, ta = (0,6717, 0,9625, 0,7247). 
Teorema de Guldin 
Sea D un dominio en el plano OXZ, que para mayor facilidad 


supondremos que es una región elemental de tipo I (el resultado 
es válido para cualquier dominio elemental): 


D=((x,y)]a<x<b, 9 (x)<2<0,(x)). 


Si es A el área de D y (x,y) su centro de gravedad, sabemos 


que 


Añ= [ideas 


Supongamos que D gira alrededor de OZ; 
engendrará un recinto 2 que será, pasando a 
coordenadas cilíndricas, el transformado de 


(para mayor claridad, en la figura hemos hecho un giro incompleto). El volumen de $2 
será 


Vs io ["raoazdr=2m [| [io rdzdr=20%A. 


Es decir, hemos demostrado el siguiente 


Teorema (Guldin). El volumen del sólido engendrado por el giro de un dominio 
elemental plano D alrededor de un eje situado en su plano y que no lo corta, es igual al 
producto del área de D por la longitud de la circunferencia descrita por su centro de 
gravedad. 


Como aplicación resulta inmediato, por ejemplo, el cálculo del volumen del toro T, 
engendrado por un círculo de radio r al girar alrededor de un eje situado a una distancia 
R de su centro (R >r): 


V(T)=xr? 27R=2x*Rr”, 


como puede comprobarse fácilmente mediante una integral simple. 
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Momentos de inercia 


El momento de inercia (o momento de segundo orden) de un sólido £2 respecto de un 
eje, el OZ por ejemplo, vendrá dado, si es mx, y, 2) su densidad, por 


q. 2 2 
[, = LG + y )m(x, y,2) dV , 
y permutando las variables obtendremos las expresiones para /, el,. 


Ejemplo 


Momento de inercia de un paralelepípedo rectángulo, de dimensiones axbxXc y 
densidad constante m, respecto del eje de simetría paralelo al lado de longitud a. 


fu ml P l (y? + 2?) dx dz dy = mal y lÑ dy = mal (o +5) dy = 
E ES > 31 =$ 12 


[22.221 a(eb' +0) — (b?+c?)M 
=ma —+— | ==, 


=S m= 
3" 1]. 12 12 
2 


Sea e una recta y G el centro de gravedad de 2. Tomemos como eje OZ una recta 
paralela a e que pase por G, y sean entonces 


x=a4 
pb 
las ecuaciones de e. Como 
2 2 2 2 2 
(x-a) +(y-b) =x?* + y? -2ax-2by+d*, 
siendo d? = a? +b? el cuadrado de la distancia entre ambos ejes, al integrar resulta 


1,=1,+M2?, 


donde M es la masa de £2, es decir, un teorema análogo al obtenido en el plano. 
Recordemos su enunciado: 


Teorema. El momento de inercia respecto de un eje cualquiera es igual al momento 
respecto del eje paralelo que pasa por el centro de gravedad más el producto de la masa 
por el cuadrado de la distancia entre ambos ejes. 
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Ejemplo 


Calculemos el momento de inercia respecto del eje OZ de la esfera vista en el ejemplo 2 
de la página 34, pero considerando ahora las longitudes en metros y la densidad en 


Tm 1,3 + Operando en coordenadas esféricas es x? + y? = p? sen? y, con lo que 
a pr px 
L.= f f f p” sen” p(2a— pcos p)p” senp dO dy dp = 
S 2 f (2ap* sen* p— p*sen* pcos p) dp do = 


= 27 ¡a (2ap*(sen p—cos? p sen p) - p” sen* pcos p) dp dp = 


3 TN 
= 27), 2 cos Qp+ e - p ce dp = 


p=0 


6 
| y _lóxa 
3 


0 
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